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An  efficient  numerical  method  based  on  the  application  of 
Galerkln's  method  In  the  spectral  domain  is  proposed  to  obtain 
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INTRODUCTION 


Monolithic  microwave  integrated  circuits  (MMIC)  based  on 
gallium  arsenide  (GaAs)  technologies  are  being  considered  seriously 
as  viable  candidates  for  satellite  communication  systems,  airborn 
radar  systems,  and  other  applications  [ 1-4  ].  The  low- loss  prop¬ 

erties  of  semi-insulating  GaAs,  combined  with  the  excellent  micro- 
wave  performance  of  metal-semiconductor  field-effect  transistors 
(MESFET's)  [5],  make  it  possible  to  design  truly  monolithic  micro- 
wave  Integrated  circuits.  By  monolithic  we  mean  an  approach  wherein 
all  passive  and  active  circuit  elements  and  interconnections  are 
formed  on  the  substrate  by  some  deposition  scheme,  such  as  epitaxy, 
ion  implantation,  sputtering,  and  evaporation. 

The  monolithic  approach  makes  it  possible  to  fabricate 
circuits  of  small  size  and  light  weight,  and  this  approach  elimi¬ 
nates  wire  bonding  between  components  in  the  circuits.  Small  size 
allows  batch  processing  of  hundreds  of  circuits  per  wafer.  The 
use  of  batch  processing  and  the  elimination  of  wire  bonding  lowers 
the  cost  of  manufacture  and  Improves  the  reliability,  reproduci¬ 
bility,  and  operating  bandwidth. 

However,  as  the  physical  dimensions  of  the  circuit  com¬ 
ponents  become  smaller,  it  becomes  more  difficult  to  trim  and 
trouble-shoot  a  working  circuit.  To  minimize  the  need  for  trimming 
and  to  simplify  trouble-shooting,  a  computer-aided  design  (CAD) 
technique  is  essential.  Another  potential  problem  that  arises  due 
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to  the  small  chip  size  Is  undesired  RF  coupling  within  the  circuit. 
To  avoid  the  coupling  problem  and  to  create  an  efficient  CAD  pro¬ 
gram  requires  a  thorough  knowledge  of  the  transmission  properties 
of  planar  transmission  lines  formed  on  semiconductor  substrates. 

(hi  a  planar  substrate,  there  are  basically  four  types  of 
transmission  lines  available,  as  shown  in  Fig.  1.  They  are 
microstrip  line  (Fig.  -1(a)),  slot  line  (Fig.  1(b)),  coplanar 
waveguide  (Fig.  1(c),  and  coplanar  strips  (Fig.  1(d)).  In 
the  past  two  decades,  these  transmission  lines  have  been  studied 
fruitfully  using  various  analytical  and  numerical  techniques  such 
as  conformal  mapping  [  6,7  ],  finite  difference  [g  J,  finite  ele¬ 
ment  [9  ],  integral  equation  [10],  spectral  domain  [11,12,13],  and 
many  other  methods  [14,15,16].  Based  on  the  considerations  of  chip 
yield  and  electrical  performance,  microstrip  lines  and  coplanar 
waveguides  are  considered  to  be  the  most  suitable  for  MMICs  [ 17] . 

One  way  of  keeping  material  costs  down  is  to  make  the 
substrate  as  thin  as  possible.  Since  thin  GaAs  substrate  is  frag¬ 
ile,  a  bottom-side  ground  plate  is  used  to  Increase  the  mechanical 
strength  of  the  coplanar  waveguide  circuits.  This  requirement 
results  in  a  conductor-backed  coplanar  waveguide  as  shown  in  Fig. 

2.  Because  of  the  thin  substrate,  the  electrical  effect  of  the 
additional  conductor  plane  backing  is  Important.  Since  no  analy¬ 
tical  results  have  been  reported  on  the  propagation  characteristics 
of  jch  a  structure,  this  is  the  first  subject  to  be  studied. 


v/n'  r  i  ’.PiO  r 


‘‘V4 


Dielectric 


OCher  structures  to  be  studied  include  slow-wave 
transmission  lines  formed  with  metal-insulator-semiconducto 
contacts  shown  in  Fig.  3.*  The  existence  of  the  slow-wav 
has  been  experimentally  observed  on  MIS  and  Schottky-barrie 
strip  lines  [18,19]  as  well  as  coplanar  waveguides  [201*  B 
of  the  reduction  in  dimensions  of  the  distributed  passive  c 
these  slow-wave  transmission  lines  could  be  advantageous  in 
design  of  MMICs.  Some  applications  employing  these  structu 
been  proposed  in  the  design  of  filters,  directional  coupler 
and  phase  shifters  [22]. 

In  the  following  chapters,  the  conductor -backed  c 
waveguide  and  the  slow-wave  transmission  lines  mentioned  ab 
be  analyzed  using  a  full  wave  analysis  in  the  spectral  doma 
Transmission  line  characteristics  such  as  propagation  const 
characteristic  Impedances  will  be  calculated  and  compared  w 
available  measurements. 

II.  FORMULATION 

The  spectral  domain  technique  was  developed  a  few 
ago  for  efficient  numerical  analysis  of  printed  transmissio 
[il. 12,13,23] •  Unlike  quasi-static  approximations,  this  me 
a  full  wave  analysis  which  can  predict  the  frequency-depend 
properties  required  in  broadband  design.  It  is  superior  to 
numerical  methods  because  of  the  efficiency  of  numerical  ca 
and  the  ease  of  formulation  [24]. 


Insulator 
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A.  Propagation  Constant 

The  cross-sectional  view  of  a  conductor-backed  coplanar 
waveguide  is  given  in  Fig.  4.  The  coordinate  system  is  defined 
in  a  manner  such  that  wave  propagations  take  place  in  the  z- 
dlrection.  The  area  of  the  cross-section  is  divided  into  two 
regions  by  the  boundary  at  y  •  d  with  each  region  containing  homog¬ 
eneous  dielectric  material. 

In  conventional  space  domain  analysis  [25],  this  structure 
may  be  analyzed  by  first  formulating  the  following  coupled  homog¬ 
eneous  integral  equations  and  then  solving  for  the  unknown  propa¬ 
gation  constant: 

/[Y  (x  -  x ' , 6) E  (x ’ )  +  Y(x  -  x',3)E  (x*)]  dx’  -  J  (x)  (  la) 

ZZ  Z  ZX  X  z 

/[Y  (x  -  x',3)E  <x')  +  Y  (x  -  x ' ,  8)  E  (x ' )  ]  dx'  -  J  (x)  (  lb) 

XZ  Z  XX  X  X 

where  E  ,  E  ,  J  and  J  are  unknown  electric  field  and  current  com- 
X  Z  X  z 

ponents  on  the  boundary  y  -  d  and  the  Green's  functions  (admittance 
functions)  Y  ,  etc.,  are  functions  of  the  propagation  constant  B. 
The  boundary  conditions  at  y  ■  d  require  that  E^  and  E^  be  zero 
except  on  the  slots  where  J  and  J  are  required  to  be  zero.  There- 
fore,  these  equations  may  be  solved  provided  that  etc*»  are 

given.  However,  for  the  inhomogeneous  structures,  the  Green's 
functions  are  not  available  in  closed  forms. 

In  the  spectral  domain  formulation,  the  convolutional 
type  of  coupled  integral  equations  (  la,b)  are  Fourier 


transformed  to  yield  the  algebraic  equations: 


?  (o,B)2f  (o)  +  ^  (a,g)2f  (a)  -  3f  (a)  (  2a) 

zz  z  zx  X  z 

Y  (a,8)£  (a)  +  ^  (o,8)if  (a)  -  J  (a)  (  2b) 

XZ  Z  XX  X  X 

where  the  quantities  with  a  tilde  ('v*)  are  Fourier  transforms  of  the 
corresponding  quantities.  The  Fourier  transform  is  defined  as 

♦(a)  r  4>(x)e:Jax  dx  (  3) 

In  addition  to  3,  the  algebraic  equations  (  2a, h)  contain  four 

other  unknowns  J  ,  Jf  ,  %  and  %  . 

Z  X  Z  X 

The  Green's  admittance  functions  are  derived  as  follows. 
First,  the  hybrid  fields  are  expressed  in  terms  of  superposition 
of  TM-to-y  and  TE-to-y  expressions  by  way  of  vector  potentials. 

The  electric  and  magnetic  vector  potentials  are  defined  in  both 
regions  (i  ■  1,2)  as 


(x.yO  *  ®y  ^(x.y) 


,-jSz 

,  i  -  1,2 

(  4a) 

,-JBs 

,  i  -  1,2 

(  4b) 

ejwt 

is  Implied,  and  harmonic 

solutions 

in  the  z  direction  are  assumed.  Beta  (3)  is  the  phase  constant 
and  a  denotes  the  unit  vector  in  the  y  direction.  Since  the  vec- 

y 

tor  potentials  satisfy  the  vector  Helmholtz  equation,  the  scalar 

Jj  6 

functions  $y^  and  should  satisfy  the  scalar  Helmholtz  equation 
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m  wji  w  ■*  -.'  ?  .  r..iv-  ^  r.  •’-  «■.  -v 


Region  I  (y  >  d): 


Yx  "  V^»2  +  *2  “  ko 

Y2  -  y«2  +  e2  -  erk^ 

where  A*1,  Ae,  and  Be  are  unknown  coefficients.  The  boundary 
condition  of  the  perfectly  conducting  plane  at  y  -  0  is  embedded  in 
this  set  of  solutions.  Also,  Is  chosen  to  be  positive  so  that 
the  radiation  condition  at  y  ■  *  is  satisfied.  Substitution  of 
these  solutions  into  the  Fourier  transforms  of  equations  (  6a, b) 

yields  the  field  expressions  in  both  regions: 


r  h  Yi  ei  -My-*) 

E*i  *  [-«*  +  J«  r  *  J  • 


-h-£l 


’ii,. 

A  e 


E.i  ■  [^ +  36  rA'l' 


’l  .e  1 


Ahj 

+  1 


’i  .hi 

i 


.  i .  ’ii.b  -^i^) 

VT,i4^lA  e 


(  Ba) 


(  8b) 


(  8c) 


(  8d) 


(  8e) 
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H 


zl 


»  r  Ahl 


-Y1<y-d) 


8f) 


Region  II  (0  <  y  <  d) 


Bx2  "  [  ^eBh  ’  j°  ^  ^  ]  8lnh  Y2y 
2 

Ey2  "[-*2  +^J®e  C°8h  V 


(  9a) 


i 

i 

i 


(  9b)  | 

i 


E 


z2 


H 


x2 


[jaBh 

[  J8B* 


e 


h 


I 

] 


sinh  y2y 


cosh  y2y 


(  9c) 


(  9d) 


(  9e) 


nz2  “  |  -JoB*  -  jB  Bh  j  cosh  Y2y 

( 

9f  ) 

Finally,  the  boundary  conditions  at  y  • 

d  have  to  be 

satisfied.  In  the  spectral  (or  Fourier  transfom) 

boundary  conditions  are: 

domain,  the 

^xl (a)  -  ^x2 (a)  -  &x(a) 

( 

10a) 

KiM  -  \zM  ■  ^(o) 

( 

10b) 

8|x2<“)  -  *xl<“)  -  Vrf 

( 

10c) 

i'J&Zt 


*  -  M 


l 


1 


V.1 


^<«)  -  s,i<«>  ■  -V»> 


(  lOd) 


*V, 

where  J  and  J  are  Fourier  transforms  of  unknown  current  components 

2  X 


J  and  J  on  the  conducting  region  at  y  ■  d.  These  conditions  allow 

2  X 


ll  6  tl  6 

us  to  obtain  the  expressions  for  coefficients  A  ,  A  ,  B  and  B  in 


•v,  <\, 

terms  of  J  and  J  : 

z  x 


A*1  *  B*1  sinh  “  j  - - - 

2  Y1  Y2 

~ —  +  ~z~  coth  y-d 

Z1  z2 


•u  ^ 
aj  -  $J 
z  x 


2  2 
a  +  B 


(  Ha) 


Y1  e  Y2  e  1 

-r=-  A  -  -  —  B  sinh  y-d  -  -j  - - x - 

yi  y2  yi  72 

—  +  —  coth  Yod 

T1  y2  2 


BJ  +  oJ 

Z  X 


2  2 
o  +  B 


(  Hb) 


By  substituting  equations  (  lla,b)  into  (  8a)  and  (  8c) ,  we 


obtain  the  algebraic  equations  (  2a, b)  with  the  Green's  admit¬ 


tance  functions  $  ,  etc.,  expressed  in  closed  forms  as  follows: 

22 


zz  2 
a 


i2  /Y1  y2  \  B2  /y  1  y2  ) 

^—jl^  +  ^coth  Y2d  +-f> +  r  coth  Y2dj 
+  yzj  z2  ^  /  o  +  b  \Yi  y2  j 


(  12a) 


Y  -  Y 
zx  xz 


■  tM?  *  ?  coth  coth  yj  < 

a2  +  S2  VI  Y2  2  *!  *2  / 


% 

Y  - 

XX 


-s—i  +  £  coth  Y2dJ  +  +  h 

*  H'l  *,  2/  o2  ♦  (2  \Yi  r2 


coth  y2<1 


(  12c) 


- - - - _____ - - - 
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This  process  of  formulating  the  Green's  functions  is 


straightforward  but  rather  lengthy,  especially  when  more  dielectric 
layers  are  used.  Therefore,  a  simpler  method  of  formulation  pro¬ 
posed  by  Itoh  [24]  will  be  discussed.  This  new  method  is  fast  and 
easy  while  it  gives  more  physical  Insight  of  the  spectral  domain 
method . 

From  the  form  of  the  Inverse  Fourier  transform 

♦  (x ,  y)  e“^Z  -  $  (a,y)  e”^0*  +  ^  do  (  13) 

y2if 

we  recognize  that  the  field  components  in  the  space  domain  are  the 
superposition  of  nonuniform  plane  waves  (spectral  waves)  propa¬ 
gating  in  the  direction  of  6  from  the  z-axls,  where 
0  »  cos-1  (8/^0^ 4  B2)  .  For  convenience,  we  define  in  Fig.  5 
a  u,v-coordinate  system  for  each  spectral  wave,  which  has  a  rota¬ 
tional  relation  with  the  x,z-coordinate  system  defined  by 

u  -  x  cosfl  -  z  sin6  (  14a) 

v  «  x  sin0  +  z  co80  (  14b) 

For  each  6  (i.e.,  for  each  a),  the  nonuniform  plane  wave  may  be 
decomposed  into  TW-to-y  ,  %v,  &u>  and  TE-to-y  (Sf^,  Eu,  Sfy)  sur¬ 
face  waves  for  which  homogeneous  boundary  conditions  apply.  Since 
the  spectral  current  is  due  only  to  the  discontinuity  of  the 
tangential  ifu  component  in  the  TM  fields,  and  3fu  is  due  only  to  TE 
fields,  they  can  be  dealt  with  Independently. 


wrnmmwmmMmmmmmmMm 
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For  surface  wave  modes,  the  characteristic  equations  may 
be  obtained  from  the  transverse  resonance  conditions  of  the  equiva¬ 
lent  transmission  line  circuits.  Fig.  5  shows  the  equivalent 
circuits  for  both  TM  and  TE  cases.  The  characteristic  admittances 
in  each  region  are 


TMi 


% 


ui 


K 


Vi 


i  -  1,2 


15a) 


Y 


TEi 


i  -  1,2 


15b) 


where  y^  «  ^a2  +  B2  -  is  the  propagation  constant  in  the  y 

direction  in  the  i-th  region.  All  the  boundary  conditions  for  the 

TE  and  TM  waves  are  readily  incorporated  in  the  equivalent  circuits. 

For  Instance,  the  conducting  plane  at  y  *  0  is  represented  by  a 

short  circuit.  The  electric  fields  E  and  if  are  continuous  at 

v  u 

y  -  d  and  are  related  to  the  currents  via 


Jv(o)  -  $e(a,B)  £v(a)  (  16a) 

Ju(o)  -  *h(a,B)  Su(a)  (  16b) 


Y  and  Y  are  the  input  admittances  looking  into  the  equivalent  cir¬ 
cuits  from  the  current  sources  at  y  •  d,  and  are  given  by 

$*(<*, 8)  -  Y®  +  Y®  (  17a) 


$h(a,B)  -  Y J  +  Y1 


(  17b) 
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where  and  are  the  input  admittances  looking  upward  and  down- 

.h 


ward,  respectively,  at  y  =  d  in  the  TM  equivalent  circuit,  and  Y 


and  Y^  are  input  admittances  in  the  TE  circuit: 


TMl 


*  Y2  *  YTM2  coth  Y2d 


(  18a) 


TE1 


»  Y2  "  YTE2  COth  Y2d 


(  18b) 


The  final  step  is  to  derive  the  relation  (  2a, b)  in  the 

x,z-coordinate  system  from  equations  (  16a, b)  via  a  coordinate 

transformation.  Because  the  transformation  is  a  simple  coordinate 

rotation  (  14) ,  E  and  E  (J  and  J  )  are  linear  combinations  of 

X  Z  X  z 

E  and  i!  (J  and  J  ).  After  the  transformation,  the  admittance 


U  V  U  V 

matrix  elements  in  (  2)  are  found  to  be 


2  %h  2 

Y  (0,6)  -  sin  e  Y  +  cos  6  Y 
zz 


(  19a) 


^  (ct,B)  ■  ^  (a, 8)  -  sinQ  cos9  ($6  -«•  Y*1) 

zx 


xz 


Y  (a, 8) 

XX  ’ 


cos2e  Yh  +  sin2e  Y* 


(  19b) 

(  19c) 


where 


sine  - 


v6*~ 


CO80 


8 


+  8*  Va2  +  82 

It  is  easily  shown  that  this  equation  is  identical  to  equation 

(  12). 


•  *  «  •  * /  »  '  «  *  •  _*•  V  V  ,4  ,4  V  ,*  *  ,•  .4  V  «•  *  - 

1  .  '  ‘  H  '4  »  •  -»  •*  *  *  ■»  s.  •  *.  *  '  *  ■  V  *  '  »  *  • 
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Now  that  we  have  derived  the  equations  correspond 

the  integral  equations  (  la,b),  with  closed  form  express! 

the  Green's  functions  in  the  spectral  domain,  the  solution 

phase  constant  B  in  the  z-direction  is  to  be  calculated.  1 

end,  Galerkin's  method  is  applied  by  expanding  the  electric 

components  G  and  G  in  terms  of  known  basis  functions  G 
r  x  z  xn 

as  follows: 


where  c  and  d  are  unknown  coefficients.  The  basis  functJ 
n  n 

and  G  must  be  chosen  such  that  their  inverse  Fourier  trar 
zn 

are  nonzero  only  in  the  slots.  After  substituting  (11-20) 

(11-2)  and  taking  the  inner  products  of  the  resultant  equat 

with  the  basis  functions  If  and  ll  ,  we  obtain  the  linear 

xn  zn’ 

* 

taneous  equations: 


l  c 

N 

+  l 

K(l,2) 

d  *0  ,  m  *  0,1. .  .M 

Ln  mn  n 

n*0 

n*l 

mn 

n  *  ’ 

y  k(2-»  c 

N 

+  l 

k(2.2) 

d  *  0  ,  m  ■  1,2. . .N 

'  _  mn  n 

n«0 

n-1 

mn 

n  *  ’ 

where 

k(1,1)  -  r  fc 

(a)  $ 

(a.B) 

£  (a)  da 

mn  *  zm  zz 

zn 

TF*.T 


'■*  -*  'A  A'A'A  'A1.1- " 


Kd,2)  _ 
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eta 
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2f 
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2! 
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l 

da 

mn 
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zm 

zx 

xn 
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(a) 

'V 

Y 

(a, 3) 

E 

(a) 

da 

mn 

—00 

xm 

xz 

zn 

k(2’2>  - 

mn 

r 

E 

xm 

(a) 

<v. 

Y 

XX 

(a, 3) 

<\< 

E 

xn 

(a) 

da 

Notice  that  the  right  hand  sides  of  (  2a, b)  are  eliminated  by 

this  procedure.  This  can  be  verified  by  using  Parseval's  theorem 
that 


r  ^  (a)  J  (a)  da  -  J*  E  (x)  J  (x)  -  0  (  22) 

zm  z  ‘  zm  z 


The  above  relation  is  true  since  E  (x)  [the  inverse  transform  of 

zm 

E^Ca)]  and  J2<x)  are  nonzero  in  the  complementary  regions  of  x  at 
y  »  d. 

Finally,  the  propagation  constant  3  is  obtained  by 
solving  the  system  of  simultaneous  equations,  (  21).  Since  this 

set  of  equations  is  homogeneous,  non-trivial  solutions  can  be 
obtained  only  when  the  determinant  of  its  coefficient  matrix  is 
equal  to  zero,  i.e.,  the  matrix  is  singular.  This  results  in  a 
characteristic  equation,  from  which  3  is  obtained. 

B.  Characteristic  Impedance 

Due  to  the  non-TEM  nature  of  the  conductor-backed 
coplanar  waveguide,  the  definition  of  characteristic  impedance  is 
not  unique.  One  possible  choice  is  to  define  it  as 


21 


Z 

c 


avg 

where  V  Is  the  slot  voltage  and  P 

o  avg 

flow  on  the  line,  which  Is  given  by 


(  23) 

Is  the  time-average  power 


Pavg  "  Re  /r  E  x  5*  •  az  dxdy  (  24) 

Since  the  limits  of  Integration  are  infinite,  Paseval's  theorem  may 
be  used  to  transform  the  expression  into  the  spectral  domain  as 


Pavg  "  Re  /r  t^x(a,y)  &*(a,y)  -  Ey(o,y)  H*(a,y)]  dxdy  (  25) 

Expressions  for  the  electric  and  magnetic  fields  in  both  regions 
are  given  by  equations  (  3)  and  (  9) .  Since  the  y-dependence 

of  these  expressions  is  simple,  the  integration  with  respect  to  y 
may  be  performed  analytically.  (  25)  then  becomes  a  single  inte¬ 

gral  of  the  form 


P  ■  R  /  g(a)  da 
avg  e  J  ° 

where 


g(a) 


-  j 


1 

2 


+  j 


i  — —  (r  ~  (“2 

sinh  Y2d  'T2  z2f  ' 


C2 


) 


(  26) 
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Note  that  the  calculation  of  these  expressions  requires  the  value 
of  if  and  il  he  known,  l.e.,  the  expansion  coefficients  c  and  d 
be  obtained  by  solving  the  equations  (  21a, b)  for  a  known  B. 

Finally,  the  voltage  VQ  is  computed.  This  Involves  sim¬ 
ply  the  Integration  of  the  assumed  electric  field  distribution 
across  the  slot  and  can  be  done  analytically.  However,  the  evalua¬ 
tion  of  (  26)  must  be  done  on  a  digital  computer. 


C.  Choice  of  Basis  Functions 

Any  kind  of  basis  function  may  be  used  as  long  as  it  is 
nonzero  only  in  the  slot  regions.  However,  due  to  the  variational 
nature  of  the  approach,  the  efficiency  and  accuracy  of  this  method 
depend  greatly  on  the  choice  of  basis  functions.  In  this  study, 
the  basis  functions  are  selected  in  accordance  with  the  following 
criteria: 


1)  For  rapid  convergence  of  the  solution  the 

functions  should  satisfy  the  edge  condition 

[  271  which  requires  that  E  (x)  behaves  like 

z 


v\y.>  .y.'/v-vv  o.-v  %•  v ■.  v-V-'V  -  .•  jk.’-'.-i 


vv 
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1  /9 

|x  -  xg|  near  the  edge  xe  of  a  strip, 
whereas  E^(x)  Approaches  xg  with  the  sin¬ 
gularity  i  x  -  xer1/2. 

2)  The  set  of  functions  E  . ,  E  .  should  be 

xl  zi 

complete  to  enable  approximation  of  the 
exact  solution  to  any  degree  desired  sim¬ 
ply  by  increasing  the  number  of  terms  of 
the  expansion.  In  this  way  the  numerical 
solutions  can  be  easily  checked  for  their 
convergence. 

3)  All  the  physical  insight  available  should 
be  incorporated  into  the  choice  of  expan¬ 
sion  functions  so  that  the  combination  of 
them  can  closely  represent  a  modal  field 
distribution  and  the  matrix  size  can  be 
held  small  for  the  given  required  accuracy. 

4)  It  is  desirable  that  the  transforms  if  . 

xi 

E^  are  available  in  a  fairly  simple  analy¬ 
tical  form. 

With  the  above  considerations  in  mind,  the  following  set  of  func¬ 
tions  are  employed  128] 


SI 


Note  that  these  functions  are  defined  only  over  one  slot.  For  even 
mode  solutions,  complete  figures  for  the  first  three  basis  functions 
are  shown  in  Fig.  7.  The  Fourier  transforms  of  the  entire  set 


(o) 

xn 


j  sinbo  [  JQ(|wa  +  y-j)  +  J0(|«  -  y|)] 


;  for  n  -  0,2,4,...  (  28a) 


-i  /f  cosba  |  J0(|wa  +  y|)  -  J0(|wa  -  yl>] 


;  for  n  -  1,3,5, . . . 


v  v  *.»  v  v  v  V >Vv' 
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cosba  | JQ(|wa  +  ^j-i>  +  Jo(|wa  -  ^yi)J 

;  for  n  ■  1,3,5,...  (  28b) 

sinba  |jo(|w»  +  ^|)  -  Jo(|wa  -  -y-|)| 

;  for  n  -  2,4,6, . . . 

where  J  denotes  the  zero-order  Bessel  function  of  the  first  kind, 
o 

D.  MIS  Slow-wave  Structures 

For  the  MIS  structures,  the  analysis  is  more  involved  due 
to  the  presence  of  a  lossy  semiconductor  layer  with  a  relatively 
high  conductivity.  The  phase  constant  and  the  bulk  attenuation  can 
no  longer  be  calculated  by  means  of  a  perturbation  theory  as  for 
the  low-loss  case.  The  losses  must  be  taken  into  account  during 

the  formulation.  This  is  done  by  employing  a  complex  permittivity 

Cf  * 

e*  -  e'  -  je"  in  the  lossy  region,  where  e"  ■  —  contains  the  con¬ 
ductivity  Information  of  the  lossy  material.  The  rest  of  the  for¬ 
mulation  follows  very  closely  the  process  described  in  Section  II, 

B.  For  convenience  of  analysis  the  cross-sectional  views  of  the 
MIS  coplanar  waveguide  and  microstrip  line  are  shown  in  Fig.  3 
together  with  their  equivalent  circuits. 

* 

The  Green's  admittance  functions  for  the  coplanar  wave¬ 
guide  are  of  the  same  form  as  (  19) ,  where  the  input  admittances 

looking  into  the  equivalent  circuits  are  given  by 
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The  same  set  of  basis  functions  as  in  (  28)  are  used  in  the 

Galerkln's  procedure.  A  complex  propagation  constant,  which  is 
composed  of  the  imaginary  phase  constant  and  real  attenuation  con¬ 
stant,  is  obtained  by  a  root-seeking  computer  program. 

In  the  analysis  of  the  MIS  microstrip  line,  the  basis 
function  expansion  is  more  easily  performed  over  the  current  com¬ 
ponents  on  the  finite  strip,  rather  than  the  electric  fields  that 
extend  to  infinity.  To  do  this,  the  formulation  process  is 


modified.  Instead  of  equation  (  2),  the  following  equations  are 

obtained 


Z  (o,6)  J  (a)  +  Z  (a,8)  J  (a) 

zz  z  ZX  X 

2  (a,B)  J  (a)  +  (a,B)  J  (a) 

XZ  Z  XX  X 

where  ,  2?  and  are 

ZZ  ZX  XZ  XX 

given  as  follows: 
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(  31a) 

(  31b) 

(  31c) 


(  31d) 


(  31e ) 


Since  the  slot  E-fields  and  strip  currents  have  analogous 
properties,  the  expressions  for  and  are  used  for  Jzn  and 
J  ,  respectively.  For  fundamental  mode  studies  only  even-symmetric 


30 


J  and  odd-synnetric  J  are  used.  Their  Fourier  transforms  are 
zn  '  xz 

Jzn(a)  -  lJo(|wo  +  mr|)  +  JQ(|wa  -  nir j ) ]  ;  n  *  0,1,2,... 


(  32a) 


Jm(a)  “  [Jq(|wo  +  nir | )  -  J0<^a  -  nir|)]  ,  n  =  1,2,3,... 


(  32b) 


Ill.  RESULTS  AND  DISCUSSION 


A.  Conductor-Backed  Coplanar  Waveguide 


The  conductor-backed  coplanar  waveguide  has  been  studied 
for  a  wide  range  of  geometric  parameters.  Fig.  9  shows  typical 
results  of  the  dispersion  property,  assuming  GaAs  as  the  substrate 
(e^  ■  13.0).  Aspect  ratios  a/d  and  w/d  are  the  parameters  varied. 
Since  the  present  structure  is  a  mixture  of  a  microstrip  line  and 
a  coplanar  waveguide,  the  properties  of  either  one  become  pre¬ 
dominant  depending  on  the  structural  parameters.  As  the  slot  width 
increases  for  a  fixed  substrate  thickness,  the  characteristic 
approaches  that  of  a  microstrip  line.  This  behavior  is  shown  in 
Fig.  9a.  On  the  other  hand,  when  the  slot  width  is  fixed  and 


the  thickness  increases,  the  behavior  approaches  the  coplanar  wave¬ 
guide  case,  as  depicted  in  Fig.  9b. 

For  a  moderate  aspect  ratio  (e.g.,  a/d  “  1/3,  w/d  ■  1/3), 


the  transmission  line  becomes  less  dispersive  than  the  corresponding 
mlcrostrlp  line  with  the  same  aspect  ratio  (a/d  *  1/3) .  This  sug¬ 
gests  that  we  may  define  the  characteristic  impedance  of  the 


Coplanor 

Waveguide 


transmission  line  as  Z  *  Z  //e  ,,,  where  e  is  the  effe 

c  o  V  eff  eff 

2 

dielectric  constant  defined  as  (B/k)  and  Z  is  the  charac 

o 

impedance  when  e  *  1.  The  impedance  Z^  computed  in  this 
plotted  in  Fig.  10.  It  is  observed  that  the  character i 
impedance  decreases  as  the  slot  widths  decrease,  while  the 
constant  is  relatively  unaffected.  Therefore,  hy  adjustin 
center  strip  width  and  slot  widths,  one  can  have  independe 
of  the  phase  constant  and  the  characteristic  impedance.  F 
suhstrate,  the  usable  practical  impedance  range  is  approx! 


15-120  0. 


The  number  of  basis  functions  required  for  accur 


largely  dependent  on  the  aspect  ratios.  For  small  w/a  and 
only  Exo  is  required  since  both  the  slot  coupling  effect  a 
effect  of  the  ground-backing  are  small.  For  larger  w/a, 
needed  to  represent  the  stronger  coupling  effect  between  t 
Finally,  when  the  thickness  becomes  comparable  to  the  slot 
more  basis  functions  are  required  for  accurate  results.  I 
study,  up  to  seven  basis  functions  have  been  used  for  some 
cases.  The  computation  time  for  such  a  solution  is  about 
on  a  CDC  Dual  Cyber  170/750. 


B.  MIS  Microstrip  Lines 


An  experimental  study  on  the  properties  of  micrc 
lines  formed  on  a  Si-SiOj  system  has  been  done  by  Hasegawt 
[18].  Their  structure  is  the  one  in  Fig.  8b  with  SiO^ 


...  • —  .  •  ...» *  ....  .  . » . *  i  *  •  ^  .  *  .  '  -  *  *  ■  »  1  ^  ^  _  *  .  •  *  .  -  St  '  **.  .  *  "A.  K1  K 


IQ.  Characteristic  impedance  of  conductor-backed  coplanar 
waveguide. 


as  the  Insulator  and  SI  ( e ■  12.0)  as  the  semiconductor.  They 
observed  that,  depending  on  the  frequency  and  the  resistivity  of 
the  substrate,  three  different  types  of  fundamental  modes  could 
appear.  By  their  physical  behavior,  these  modes  were  classified 
[18, Into  the  skin-effect  mode,  the  lossy-dielectric  mode,  and 
the  slow-wave  mode.  In  particular,  the  slow-wave  mode  was  found  to 
propagate  within  the  resistivity-frequency  range  suited  to  mono¬ 
lithic  circuit  technology.  Extensive  experimental  results  of  the 
propagation  constant  and  characteristic  impedance  were  reported  for 
microstrip  lines  with  narrow  and  wide  strips  operating  in  the  slow- 
wave  mode. 

The  validity  of  this  analysis  is  verified  by  Hasegawa's 
measurements.  For  the  case  of  a  narrow  strip  (w/d  «  0.64),  the 
calculated  slow-wave  factors  and  attenuation  constants  are  plotted 
in  Fig.  11  versus  frequency.  Although  only  one  basis  function  is 
used  to  approximate  the  strip  current,  good  agreement  is  found 
between  the  analytical  and  experimental  results.  The  most  notable 
discrepancy  in  the  attenuation  constant  toward  the  lower  end  of 
frequencies  is  due  to  the  fact  that  no  conductor  loss  is  considered 
in  the  analysis.  This  is  verified  by  including  the  conductor  losses 
of  the  ground  plane  in  the  analysis.  The  results  of  doing  this  are 
presented  in  Fig.  11  as  the  dashed  lines  which  show  an  increase 
in  the  attenuation  constant  at  lower  frequencies.  The  conductor 
losses  on  the  strip  have  not  been  considered  due  to  difficulties  in 


the  formulation. 
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At  a  fixed  frequency  of  0.5  GHz  the  structure  is  studied 
by  varying  the  resistivity  of  the  semiconducting  substrate.  The 
results  are  shown  in  Fig.  12.  The  optimum  resistivity  for  low- 
loss  propagation  is  found  to  be  about  1.5  ft-mm,  same  as  the  optimum 
resistivity  predicted  using  a  parallel-plate  transmission  line 
model  118]. 

For  a  wide  strip  (w/d  -  1)  more  basis  functions  are 
required  to  obtain  good  results.  The  convergence  of  the  solutions 
has  been  tested  by  increasing  the  number  of  basis  functions  and  is 
tabulated  in  Table  1.  It  is  found  that  good  results  can  be 
obtained  by  using  five  basis  functions  (M  ■  N  *  2) .  The  slow-wave 
factors  and  the  attenuation  constants  shown  in  Fig.  13  are 
calculated  using  these  five  basis  functions.  Although  the  calcula¬ 
tions  do  not  show  close  quantitative  agreement  with  the  measurements. 


they  predict  correct  behavior.  The  discrepancies  are  caused  in 
part  hy  neglecting  the  conductor  losses  in  the  system.  Another  pos¬ 
sible  source  of  errors  may  be  due  to  dimensional  uncertainties. 

For  example,  if  we  take  the  thickness  of  the  insulator  (d)  as  0.7  pm 
instead  of  1  pm  in  one  of  the  calculations,  the  result  agrees  much 
better  with  the  measurement  as  shown  in  Fig.  13. 

The  Schottky  contact  microstrip  line  tested  hy  JSger  and 
Itabus  l  199  was  also  analyzed.  In  such  structures,  the  depletion 


'-wave  factor  and  attenuation  of  MIS  iticrostrip  lines 
resistivity  of  seaiconductor  substrate. 
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lines  pass  through  the  depletion  layer  and  the  field  at  large 
distances  from  the  depletion  layer  has  little  effect  on  the  propa¬ 
gation  characteristics.  It  should  be  noted  that  the  frequency- 
dependent  property  and  f ringing-field  effect  of  the  structure  are 
correctly  incorporated  in  the  present  method.  Such  is  not  the  case 
in  the  parallel  plate  approximation  used  in  ti8,19i  in  which  some 
correction  factor  must  be  used.  Fig.  14  shows  that  the  present 
analysis  can  accurately  describe  the  propagation  properties  experi¬ 
mentally  obtained. 

C.  MIS  Coplanar  Waveguides 

The  MIS  coplanar  structure  in  Fig.  -8a  is  also  studied 
with  this  method  and  the  results  are  compared  with  measurements 
reported  in  [20,29]  in  Fig.  15.  Good  agreement  is  shown  in  both 
the  slow-wave  factors  and  attenuation  constants.  The  coplanar  line 
was  built  on  a  GaAs  substrate  with  a  narrow  center  strip  and  wide 
slots.  Since  the  coupling  between  slots  is  strong,  a  large  number 
of  basis  functions  are  required  for  accurate  results.  The  conver¬ 
gence  of  solutions  has  been  studied  by  increasing  the  number  of 
basis  functions  used.  Typical  results  are  tabulated  in  Table  2. 
It  is  found  that  good  results  can  be  obtained  by  using  seven  basis 
functions  that  represent  axial  and  transverse  electric  fields  in 
the  slot. 
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Parameter:  Reverse  Bias 
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Slow-wave  factors  and  attenuation  constants 
of  MIS  coplanar  waveguide. 
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a(dB/mn) 

A  /A 
g 

ot(dB/mm) 

0  .019540 

.18276 

.032517 

2.8792 

1  .023043 

.060137 

.030794 

2.4416 

2  .025519 

.049147 

.031861 

1.9117 

3  .026336 

.046614 

.031633 

1.7648 

45 


Using  the  structural  parameters  given  in  Fig.  15,  a 
family  of  curves  was  generated  relating  the  propagation  constant 
of  the  MIS  coplanar  waveguide  as  a  function  of  frequency  to  the 
substrate  resistivity.  A  mode  diagram  made  from  these  curves  is 
shown  in  Fig.  '  16a.  The  three  types  of  modes  found  by  Hasegawa 
[18]  are  clearly  shown  in  this  figure.  A  typical  example  demon¬ 
strating  the  behavior  of  the  slow-wave  factor  and  attenuation  con¬ 
stant  at  0.1  GHz  with  respect  to  resistivity  is  given  in  Fig. 

16b  and  16c.  For  this  structure  the  optimum  resistivity  for  low 

_2 

loss  propagation  is  about  1.5  x  10  fi-m m. 

A  GaAs  coplanar  structure  with  Schottky  contacted  center 
strip  was  tested  at  Hughes  Aircraft  [30] .  The  measured  slow-wave 
factors  agree  reasonably  well  with  those  theoretically  predicted  by 
the  present  analysis. 

IV.  CONCLUSIONS 

An  efficient  numerical  method  has  been  proposed  to  obtain 
the  propagation  constant  and  characteristic  impedance  of  planar 
transmission  lines  used  in  MMICs.  This  method  is  based  on  the 
application  of  Galerkin's  method  in  the  spectral  domain  in  which 
the  dyadic  Green's  functions,  in  the  spectral  domain,  are  derived 
using  a  simple  and  fast  procedure.  This  procedure  is  based  on 
transverse  equivalent  transmission  lines  for  a  spectral  wave  and  on 
a  simple  coordinate  transformation.  Further,  a  complete  set  of 
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Fig.  16.  Characteristics  of  MIS  coplanar  waveguide,  (a)  Mode  chart, 
(b)  Attenuation  vs.  resistivity,  (c)  Slow-wave  factor  vs.’ 
resistivity. 


basis  functions  is  employed  so  that  the  convergence  of  sol 
can  be  tested. 

The  transmission-line  properties  of  the  conductc 
coplanar  waveguide  are  analyzed.  Since  it  is  a  mixture  of 
strip  line  and  a  coplanar  waveguide,  calculations  show  tha 
erties  of  each  become  predominant  depending  on  the  structu 
meters.  For  a  fixed  substrate  thickness,  the  character ist 
ance  and  phase  constant  may  be  varied  independently  by  sin 
adjusting  the  widths  of  the  center  strip  and  slots  in  the 
sion  line.  For  a  GaAs  substrate,  the  usable  practical  imp 
range  is  approximately  15-120  0. 

Microstrip  lines  and  coplanar  waveguides  formed 
and  Schottky  barrier  contacts  are  also  analyzed.  The  pret 
a  highly  lossy  semiconducting  substrate  is  taken  into  accc 
employing  a  complex  permittivity.  The  final  solution  is  t 
propagation  constant  that  is  composed  of  a  slow-wave  facte 
attenuation  constant.  Depending  on  the  frequency  and  the 
vity  of  the  semiconductor  substrate,  three  different  type* 
damental  modes  are  predicted.  The  optimum  resistivity  foi 
slow-wave  propagation  can  be  obtained. 
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